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Abstract 

For the N— plets of bound states in a quasi-exactly solvable (QES) toy model (sextic 
oscillator), the spectrum is known to be given as eigenvalues of an N by N matrix. 
Its determination becomes purely numerical for all the larger N > Nq — 9. We 
propose a new perturbative alternative to this construction. It is based on the fact 
that at any N, the problem turns solvable in the limit of very large angular momenta 
£ — > oo. For all the finite I we are then able to define the QES spectrum by convergent 
perturbation series. These series admit a very specific rational resummation, having 
an analytic or branched continued-fraction form at the smallest N = 4 and 5 or TV = 6 
and 7, respectively. It is remarkable that among all the asymptotically equivalent 
small expansion parameters /i ~ l/(£ + (3), one must choose an optimal one, with 
unique shift (3 = /3(N). 
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1 Introduction 



The D— dimensional and central Schrodinger bound-state problem 

[-A + V(\r\)]V(r) = £tf(f) 

finds its applications as an approximate model in nuclear physics [1] as well as in 
quantum chemistry [2] and atomic physics [3]. Last but not least, it is encountered 
as a popular exercise in textbooks on quantum mechanics [4]. In all these cases, the 
main mathematical merit of the model (1) lies in its separability, i.e., reducibility to 
the set of the ordinary (so called radial) differential equations 

— + 1 T ; + V(r) 4>(r) = £ 4>(r) , i = (D - 3)/2 + m , m = 0, 1, . . . . 
ar A r z 

(1) 

They are all easily solvable by any standard numerical technique [5] . 

Whenever we need a non-numerical solution, we must restrict the class of our 
potentials V(|r|) to a subset which is exactly solvable (ES, say, harmonic-oscillator) 
or at least quasi-exactly solvable (QES, see, e.g., the monograph [6] for a thorough 
review). Of course, the real appeal of any ES or QES potential derives from an 
agreement of its shape or spectrum with some "more realistic" numerical model 
V^ hys \\r\). In this sense, many practical applications of the ES subset are marred 
by its comparatively small size and, hence, by a very narrow variability of the avail- 
able shapes of V'^dfl) (cf., e.g., ref. [7] for their representative list). In parallel, 
the wealth of the forms of the multi-parametric wells V^^d^l) may happen to be 
more than compensated by the much more complicated tractability of their pertur- 
bations [8]. Moreover, even the very zero-order QES constructions may become quite 
complicated, whenever the size N of the required QES multiplet leaves the domain 
of N < N where the practical determination of the spectrum remains non-numerical 
[9]. Although the latter remark sounds like a paradox, one must keep in mind that 
even in the simplest QES models, the explicit construction of the larger QES mul- 
tiplets with N > N requires the more or less purely numerical diagonalization of 
asymmetric N by N matrices (see below). In purely pragmatic terms, many closed 
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QES constructions [10, 11] are in fact more difficult than an immediate numerical 
solution of the differential Schrodinger eq. (1) itself. 

In what follows we shall argue that it is possible to avoid the latter brute-force 
methods even in the QES domain where N > N . We shall recommend the use 
of the so called large— £ (or large— D [12]) expansion philosophy [13] in its new im- 
plementation which takes into account several specific features of QES models. We 
shall reduce all the inessential technicalities to a necessary minimum and restrict our 
attention just to one of the most elementary illustrative QES examples 

V [QES] (\r\) =a\r\ 2 + c\r\ 6 , c = 7 2 > (2) 

which is known to possess an elementary N— plet of QES solutions at any N [14]. 
On the background of this example we shall start our considerations in section 2 by 
showing that due to the absence of the quartic term in eq. (2), the QES solvability 
conditions are transparent and instructive (cf. subsection 2.1). In subsection 2.2 we 
show how these conditions simplify even the standard large— £ perturbation expan- 
sions of the spectrum where one uses just an asymptotic, £ ^> 1 simplification of the 
differential Schrodinger equation (1). In section 2.3 we add that in the more specific, 
QES-related setting, the algebraic construction of the solutions might be preferrable, 
at least due to its amazingly transparent form. 

The latter expectations are more than confirmed in the subsequent section 3. 
After a step-by-step analysis of the QES secular equations for the energies in depen- 
dence on the growing dimension N, we are able to reveal their general polynomial 
structure and perturbation form. The first subsection 3.1 shows that the models 
with N < 3 are all fully solvable via an elementary re-scaling of the matrix secular 
equation. Formally, our presently conjectured "optimally shifted" large— £ series re- 
main trivial in this case. The zero-order result is exact and its secular polynomial 
does not contain any perturbation. All the higher-order corrections simply vanish. 
Still, the structure of the N = 3 secular equation is already rich enough to illustrate, 
why the "optimal" value of the shift of £ 3> 1 is so exceptional, and how a return to 
its "non-optimal" values would re-introduce non-vanishing perturbation corrections. 
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In subsection 3.2, some of these observations are re-confirmed at N = A. Ele- 
mentary analysis offers the explanation why our "optimally shifted" large— I power 
series are convergent. The determination of the explicit value of the radius of conver- 
gence fj, m ax reveals that the circle of convergence covers the whole domain of physical 
interest, i.e., formally, all the positive spatial dimensions D > 0. 

Subsection 3.3 uses the slightly modified N = 5 example and shows that a Pade- 
like re-summation of our power series exists. The "re-summed" energies E (as well 
as their squares Q = E 2 ) acquire an analytic continued-fraction form. 

In subsection 3.4, the "first nontrivial" example is shown to emerge at N = 6. 
We sample there the energies in their "generalized" or "branched" continued-fraction 
form as well as in their more standard power-series representation. A few questions 
concerning their convergence are clarified. 

In subsection 3.5, an extrapolation of all the preceding constructions to any size N 
of the QES multiplets is described and discussed. For quadruplets of N = AK, AK + 
1, Ak+2 and AK+3, our main result is formulated as an iterative K— term formula for 
the systematic construction of the energies in a generalized continued-fraction form. 
Simultaneously, the Taylor-expansion technique is shown to generate the standard 
power series with the "optimally shifted" large— I perturbation structure. The last 
illustrations are added showing the smoothness of transition from the "last solvable" 
N < N = 9 to "the first unsolvable" N = 10. 

In section 4 a thorough re-interpretation of our optimally shifted large— t per- 
turbation series is given in the more standard language of the textbook Rayleigh- 
Schrodinger perturbation theory. Once more we summarize our approach in the 
last section 5, as a source of new types of perturbation expansions tailored for the 
computation of the QES spectra at any multiplet size N. 
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2 Conditions of quasi-exact solvability 



2.1 Toy model: sextic oscillator 

Sextic oscillator (2) may support elementary bound states 



JV-1 



1 



provided only that we choose one of the couplings in consistent manner [14], 



a = a (N) = -7 (AN + 2i + 1) 



(3) 



(4) 



The latter condition increases a phenomenological appeal of our example by assigning 
to potential (2) the manifestly non-perturbative double-well shape. Many compu- 
tational difficulties of "realistic" calculations may be mimicked, even within such a 
framework, whenever one chooses a larger degree N in wave function (3). Integer 
N measures the "dimension" of our problem since, under the constraint (4) and af- 
ter the insertion of (3), our differential Schrodinger equation (1) is transformed into 
equivalent matrix problem 
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(5) 



y h-N-i J 



The matrix of this system is asymmetric, 



A n = -4 7 (iV-n), C n = -2(n+ 1) (2n + 2£ + 3), n = 0, 1, . . . (6) 

so that up to the first few lowest dimensions N, equation (5) need not be easy to 
solve at all. The advantages gained by the polynomiality of ip may be quickly lost 
with the growth of their degree N. The merits of the "exceptional" QES levels 
seem to disappear in the domain of the larger N 1. However, these states may 
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"remember" and share some features of the completely solvable harmonic oscillator. 
In this sense, our oversimplified toy example may serve as a guide towards a future 
better understanding of the partially solvable Schrodinger equations with more free 
parameters [11]. 

2.2 Large— I domain in coordinate representation 

For our particular QES model of section 2.1, the radial Schrodinger eq. (1), 

^(r)=£^(r), (7) 



d _ 7(4Ar + 2£+1)r 2 + 7 2 r 6 



dr 2 r 

may be treated by the so called 1/Af (i.e., in the present notation, l/£) expansions 
along the lines recommended in the recent review [15]. This approach visualizes the 
Hamiltonian in eq. (7) as composed of an "effective" kinetic energy in one dimension 
and of an "effective" one-dimensional potential well V e ff(r) with a steep and repulsive 
peak near the origin and also with a confining growth at large r. Its minimum is 
unique and may be localized at the point 

r = arc = +*(*> )''', (8) 

A(£) = VlQN 2 + 16NI + 8N + 1Q£ 2 + 161+ 1. 
At an illustrative strength 7 = 1 and for the simplest N = we have 

V eff (r) = i^— -2£r 2 + r^ [l + 0(l/£)} = 16 £ [r - R(£)} 2 + 16 £ 3/4 [r - R(£)f + . . . . 

As long as the position of the minimum grows with £ (or D) to its asymptotically 
leading-order value R(£) = {£/^) l / 4 {\ + 0(1/ £) (note that we switched to 7 = 1 
temporarily), we may abbreviate r — R(£) = £/(2£ 1 / 4 ) and replace our Schrodinger 
equation (7) by its £ ^> 1 asymptotically equivalent anharmonic-oscillator form 



de 2yf£ 16£ U 3/2 



0(0 = ^0(0 0) 
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In terms of perturbation theory, the leading-order energies coincide with a re-scaled 
equidistant spectrum of harmonic oscillator, 

£ = AVI(2n + 1) [1 + 0(1/ VI)], n = 0,l,... (10) 

All the higher-order corrections may be computed in systematic manner, leading to 
an asymptotic, divergent [15] power series in A = 0(1/ Vt)- In our present paper, 
we are going to replace them by convergent series in the powers of A 4 . 

2.3 Large— £ domain in matrix representation 

Partially, our present paper has been motivated by all the realistic models (1) where 
the dimension D is large enough [like, e.g., in ref. [1] where D = O(10 3 )]. In the 
preceding subsection we also saw that there might exist correlations between the 
changes of D (or £) and of the conditions of (quasi-) exact solvability. Let us now 
pay more attention to this phenomenon which appeared as a source of inspiration 
in several studies of QES problems [16] where it has been noticed that, at a fixed 
matrix dimension N, the solutions are getting simpler whenever the spatial dimension 
D grows to infinity. 

In an introductory step we imagine that the QES constraint (4) merely replaces 
the definition (6) of the lower diagonal in eqs. (5) by the shorter formula A n = 
47 (n — N). An uncomfortable asymmetry of our linear algebraic problem Q(£)h = 
(5) is still there but it may be weakened by the shift of the large integer 2£ = 
D + 2m — 3 (in the m— th partial wave) to another (and also large) auxiliary integer 
G = D + 2m + N — 2. This shift represents Q(£) as a slightly more balanced matrix 
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dominated by its upper diagonal, 

( £ 2(G + 2-N) 
4(iV-l) 7 £ 4(G + 4-iV) 



87 



V 



£ 2(N - 1)(G - 2 + N) 
4 7 £ 



We re-scale it by its pre-multiplication by the diagonal matrix p with elements pP 
where p - 



'G/(2 r y). The parallel renormalization of the Taylor coefficients hj 
Pj = [G/(2^)]^ 2 hj enables us to re-scale the energy, 



£ = 2J2-fGE, 



(11) 



re-expressing QES spectrum in terms of eigenvalues of asymmetric two-diagonal ma- 
trix 

fM 

/(iv-i)(0) 



H{p) = 

\ 

where we abbreviated 

f n (p) = n - (N - 2n)/j,n, 



/2M 
/ 2 (0) 





A(0) 



f(N-l)(») 



1 



(12) 



» = 77 = 



(13) 



G 2m + N + D-2 

This matrix will now be understood as our QES "Hamiltonian" which is to be diag- 
onalized perturbatively. In a way guided by a few experiments made at the first few 
lowest dimensions N, we shall reveal and describe a new and efficient recipe for such 
a construction in the next section. 
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3 Perturbation expansions in /j 2 



3.1 Optimal shift: matrix dimension N = 3 



In our notation where the spatial dimension D and all the angular momenta I = 
(D — 3)/2 + m with m — 0, 1, . . . are "very" large, D 3> 1, the integer N denotes 
just the matrix dimension in eq. (5) and its numerical value may be arbitrary. Thus, 
our perturbation expansions will use the smallness of // = 1/G or, in a more general 
setting, of // = 1/[D + 2m + S(N)] where we have a full freedom in the choice of an 
iV-dependent shift 5 = 5(N). 

Let us select 5(N) = 1 + e at iV = 3, and reduce the diagonalization of our 
"Hamiltonian" (12) to the elementary secular equation 



det 



E l-fi(l + e) 

2 E 2 + 2ji(l-e) 

1 E 



0. 



Although its explicit polynomial form E 3 — AE + AEfie = 0is virtually trivial, it 
may be still simplified and made perturbation independent, provided only that we 
choose the above shift 5 with the special and unique e — 0. After this "optimal" 
choice, all the three zero-order energy roots Ej e {—2,0,2} remain safely real and 
become manifestly independent of the spatial dimension D. As a consequence, all 
the higher-order corrections vanish at e — 0, and the only D— dependence of the 
energies S remains encoded in the above-mentioned scaling rule (11). 



8 



3.2 Continued fractions: matrix dimension N = 4 



det 







In the next case at N = 4 one reveals that with another "optimal" 5(N) = 2, an 
explicit linear //— dependence disappears again from the secular equation 

E 1-2/i 
3 E 2 

2 E 3(1 + 2//) 
1 E 

Only the quadratic term persists in the secular polynomial E 4 — 10 E 2 + 9 — 36 /j, 2 . 
Thus, the choice of N = 4 should be understood as "the first nontrivial" model. Not 
nontrivial enough of course: the quadruplet of the secular roots keeps its closed form, 








Ej e {±y 5 ± 2 a/4 + 9 /i 2 }. Circumventing the standard and tedious perturbation 
recipes, all these four roots may be Taylor-expanded giving directly all the four 
perturbation expansions in the powers of fi 2 , 



33 



309 



E 1 = -E 4 = 3 + -if - —fT + —fi 



64' 



9 



512' 
2187 



14133 
" 16384 
137781 



// 



-fi - 



-fi - 



179643 
131072^ 
3601989 



10 + O(// 2 ), (14) 
//° + 0(/, 12 ) . (15) 



E 2 = -E 3 = 1 - -ii 2 IT - 

6 4 P 64 P 512 r 16384 r 131072 

Due to the analyticity of these expressions in the complexified variable fi, all of these 
four series have, obviously, the same and non- vanishing circle of convergence with the 
radius fi ma x = 1/2 determined by the nearest branch point in the complex plane of 
fi. Once we recollect that for any m > and D > we always have 1//J, > 5(4) = 2, 
we may conclude that for all our four shifted-large— I perturbation series at iV = 4, 
the circle of convergence is amply sufficient to cover the whole domain of physical 
interest. 

The overall square-root form of the "energies" E might have been taken into 
account in this context. At least, one may expect a simplification of the formulae 
when the squares of the energies Q = E 2 are taken into consideration. For example, 
in place of the latter two series (15), the "unified" power-series formula reads 



+ O (//) . 



(16) 
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Its use improves the rate of convergence because we removed one of the branch 
points, arriving at the new value fi [ ZZ roved] = 2/3 for the radius of convergence. In a 
broader context, it is therefore not surprizing that the underlying idea of expanding 
a suitable function of E (in our case, a square) rather than this observable quantity 
itself did find non-trivial applications elsewhere [17]. 

The transition from E to its square ft inspires also a slightly counterintuitive 
re-arrangement of our secular polynomial in an incompletely factorized form where 
the small corrections are manifestly separated, 

(ft -9) (ft-l) = 36/x 2 . 

This defines an explicit analytic continued-fraction re-arrangement of the complete, 
"perturbed" ft(/x) = ft(0) + /i 2 cj(/i). For example, once we select ft(0) = E\^ = 9, 
we arrive at the easily verified identity 

ft( / i 2 ) = 5 + 2 v /4 + 9/i2 = 9 + -|^ z . (17) 

° + s+... 

The backward conversion of the analytic continued fraction on the right-hand side 
to the plain perturbative power series in /i 2 is routine [18] and almost as trivial as 
its above-mentioned Taylor-series derivation from the left-hand-side square root. 

3.3 Convergence of iterations: matrix dimension N = 5 

Our experience gained from the preceding, completely and easily solvable examples 
applies to all of the higher iV's in a more or less straightforward manner. Firstly, by 
extrapolation we conjecture (and subsequently verify) that the "optimally shifted" 
small parameter \i is always defined by eq. (13). Secondly, we pick up the next 
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N = 5 and re-write the secular equation 



det 



E 1 — 3/2 

4 E 2(1 — yu) 
3 E 3(1 + ji) 













2 




E 4(1 + 3 n) 
1 £ 



= 



in the polynomial form E 5 — 20 E 3 + 64 E — 288 E fj 2 = showing that we may factor 
the zero root E = E 3 = out. Later on we shall see that such a root exists (and will 
be dropped) at all the odd N. The fact that one of our roots is identically vanishing 
and may be factored out, strengthens the parallels between N = 4 and N = 5. In 
the latter case we get a quadruplet of the nontrivial roots, 

75957 



I , 57 030 

E 1 = -E 5 = ^/l0 + 6 v /l + 8/i 2 = 4 + 3/i 2 -— /i 4 + " 



32 



/! - 



E 2 = -E 4 = ^10 - 6 ^1 + 8/i 2 = 2-6/i 2 +3 /U 4 -39/i 6 + ^ 8 -^/ +O (/i 12 ) . 

In full parallel with iV = 4 it is easy to show that [x ma x — 1/(2 v^2)- This time, the do- 
main of the perturbative convergence covers more than we really need. On its bound- 
ary we set m = and D = 1 and get \i = fi max phys = 1/4 < 1/(2 \/2) = Umax math- 
In comparison with N = 4 & tendency is observed towards and improvement of the 
overall rate of convergence, which are good news in the light of our intention to move 
towards any N in principle. 

Let us once more recollect the existence of the N = 4 continued fractions for fi's. 
It is related to the partial factorization of the secular equation which may be easily 
transferred to iV = 5, 

(fi-16)(fi-4) = 288/i 2 . (18) 

Re-scaling the measure of perturbation /z 2 = A/72 and picking up the zero-order 
root = 4 for definiteness, the corrections in Q(li 2 ) — 4 (1 — Z) may be defined 



512 
3693 



^ + O (/°) , 
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by iterations of eq. (18), Z[ neu) ] = A/(3 + Z[ i^), and they coincide again with an 
analytic continued fraction, 

Z = (19) 

The proof of its convergence is easy when performed in the spirit of the fixed- 
point method of ref. [19]. Indeed, the mapping — > Z[ new } has just two 
fixed points, the values of which are known and given by the quadratic equation 
Z|ppj = Z[ i d ]) = Z[ new ]. In the next step one imagines that in the vicinity of these 
two points, the obvious sufficient condition for the convergence/divergence of intera- 
tions of our mapping (i.e., of our continued fraction with constant coefficients) is that 
the derivative Y(z) = Z| nOT j (of the dependent variable Z = Z[ new ] with respect to 
the independent variable z = Z\oid\)) is smaller /bigger than one in its absolute value, 
respectively [19, 20]. We have all the necessary quantities at our disposal so that 
the criterion is verified by trivial insertions. In perturbative regime of the smallest 
A 1 it is sufficient to specify the leading-order estimates for the small positive 

Z\fp] — A/12 + . . . and for the much larger and negative -Z^p] = —3 — A/3 + We 

may conclude that the large fixed point is unstable since = — 9/A + . . . <C — 1 
while the small fixed point is stable since Y^ = —A/9 + ... E (—1, 1). 

We described this proof of convergence of continued fractions in full detail because 
its idea is easily transferred to all the higher dimensions N. In the language of fixed 
points (or accumulation points), one can also better understand the transition to the 
power-series perturbation expansions of our squared energies Q(fi 2 ). They may be 
re-constructed directly from the recurrences [20], without any recourse to the closed 
square-root formulae. 

3.4 N = 6 and the branching continued fractions 

At N = 6 the extraction of the roots from the secular equation 

-225 + 259 E 2 - 35 E 4 + E 6 + 3600 /i 2 - 1296 fi 2 E 2 = (20) 
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ceases to be easy. Fortunately, the D — > oo asymptotic analysis of this equation 
reveals that there still exists an elementary partial factorization of this equation, 

(E + 5)(E + 3)(E + l)(E-l)(E-3)(E-5) = 144 qm 2 (3 E - 5) (3 E + 5). (21) 

Both the existence and the simplicity of this relation opens to us new horizons. 
We realize that we may succeed in an algebraic determination of the perturbation 
representations of the roots in an almost as closed form as above. What is now vital 
is that the construction does not still require any explicit knowledge of the roots and 
might remain feasible at the higher JV's. 

Once we return to eqs. (17) and (19) for guidance, we find that a specific gener- 
alization of the continued fractions exists at N = 6 since Firstly, perturbation anzatz 
fi(/z 2 ) = 1 + Ai 2 ^(/W 2 ) converts equation (21) into iterative recipe 

_ 144 (25 - 9 fl) o-i+^W ,_ni (22) 

In comparison with the continued-fraction predecessors of this nonlinear two-term 
recurrence, just a "branching" of the continued fraction occurs, 

tjW =72 f 1 — n , ^-ttt I , A; = 0,1,.... (23) 

\8-^w [oj4 j| 24- fi u [old] / 

These formulae generate the sequence of the values uj^^ = oo^lw] w hh k — 0, 1, . . . 
starting from the initial = and giving 

i2 ' ( 2 + 3 /U 2 )(2 + /i 2 )' "" 1 j 

In the limit fc — > oo (whenever it exists), one gets the function a;fer = wfcL of 
/i 2 whose form is a generalization of the standard analytic continued fraction. Its 
dependence on fi 2 > is smooth. From the negative value uj(0) = —12 in the origin 
this function decreases to a minimum cu(0. 400186) = —24.94120 and then it very 
slowly grows again. 

Whenever necessary, one can establish a number of parallels between iterations 
(23) and the current theory of continued fractions [18]. In particular, the fixed-point 
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pattern of the proofs of convergence of ref. [19] applies to functions fi(/i 2 ) at both 
N = 4 or 5 and N = 6 (and, as we shall see below, 7). Thus, for the "first nontrivial" 
N = 6 mapping — > w [n<L]; we ma y s kip the details (like, e.g., all the small— A 
analysis) and emphasize only that the proof of the convergence of our branched 
continued fraction consists of three steps now. Firstly, one demonstrates that there 
is just one fixed-point root which is compatible with the perturbation smallness of 
corrections. Secondly, one confirms that this "only acceptable" fixed point is stable 
(by showing that the derivative of the mapping at this point is sufficiently small, 
|Y(z)| < 1). Finally, one shows that the stable point is unique, which follows from 
the observation that > 1 at the other two fixed points. 

For giving to the reader a rough estimate of the rate of convergence of the it- 
erations, let us pick up a sample value of the parameter A = /i 2 = 1/10, giving 
the fairly small Y(z) w 0.25949 at the accumulation point z —2.30666, while the 
quick divergence of the iterations at the other two fixed points results from the large 
magnitude of the derivatives F(5.26698) « 3.0556 and F(37.03969) « -5.5675 there. 

For many practical purposes, it is desirable to convert our branched continued 
fractions w^l = wj^L in example (22) into the more common Rayleigh-Schrodinger- 
type power series. We performed such a conversion giving, for the squared energy, 
the following perturbation series in /i 2 , 

9 , 591 fi 4215 „ 286293 in 

n = i- - 57^ 4 - — ^ + — / + -^j-^ 10 + 

3702951 12 63786951 14 3242255193 1fi 

H a 12 u 14 /i 16 - (25) 

256 ^ 1024 ^ 4096 ^ 

32707656915 1S 1182033909831 9n ^ / 99 



16384 ' 65536 
Up to the order O (/x 26 ) we verified the validity of an empirical rule that the signs are 
changing after every third order. A posteriori, one may expect an efficient numerical 
summability of this series by the standard Pade-resummation techniques. We skip 
here the really amusing possibility of a backward comparison of the resulting rational 
approximations with their available initial versions (24). 
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3.5 Expansions at an arbitrary matrix dimension N 

In a way which is outlined in Table 1 and which may be tested on the further 
non-numerical QES constructions (with dimensions N = 7 - N = 9, i.e., up to 
quartic-polynomials) as well as at any higher N, we have reached the stage where 
the overall structure of secular polynomials is clear, 



The routine analysis confirmed our expectations that at any N, the large parameters 
£ or D entering the measure of smallness /j, as defined by eq. (13) are "optimally 
shifted". This assertion is supported by the following four reasons at least. 

• All the odd powers of /x disappear from the polynomial forms of the secular 
determinants as well as from the perturbation expansions of the observable 
quantities (energies). 

• Up to the dimension N = 9 which is already fairly large, the degree of our 
"optimalized" secular polynomials does not exceed four; this means that their 
factorization may be made non-numerically. 

• At all the larger N > 9, a partial (i.e., asymptotic, D — > oo) factorization is 
still feasible non-numerically in zero order. As a consequence, all the higher- 
order C(/i 2 ) corrections may be generated in the form which, in many a respect, 
generalizes the analytic continued fraction. 

Surprizingly enough, the use of our present generalizations of continued fractions 
has led to a very transparent alternative to Cardano formulae at dimension as low 
as iV = 6. Similarly, at iV = 7 we eliminate the trivial root E = and get the rule 



V(n) = V {Q)Asy + /i 2 Pi(^) + • • • + tf K V K {®) , K = entier — . (26) 



-2304 E + 784 E 3 - 56 E b + E 7 + 40320 fx 2 E - 4320 fi 2 E 3 = 



and the parallel update of eq. (21) 




(27) 
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We need not repeat the extraction of the conclusions similar to the ones in section 3.4. 
A not too dissimilar remark may be added concerning N = 8 with 

(p - i 2 ) (p - 3 2 ) (p - 5 2 ) (p - 7 2 ) = 360 p (1225 - 682 n + 33 fi 2 ) +i5876oo p . 

(28) 

This equation illustrates both the emergence of the higher powers of yU 2 on the right- 
hand side and the parallels with the N = 9 secular equation 

(p - 2 2 ) (n - A 2 ) (p - 6 2 ) (p - 8 2 ) = 

= 288 /i 2 (22016 - 3740 D, + 99 fi 2 ) + 24385536 fi 4 

as well as an unexpectedly smooth character of transition to the "first unsolvable" 
N = 10 example 

-893025 + 1057221 E 2 - 172810 E A + 8778 E 6 - 165 E 8 + E w + 

+71442000 n 2 - 48647664 fi 2 E 2 + 3809520 fi 2 E 4 - 61776 fi 2 E 6 - 

-914457600 // 4 + 199148544 /i 4 £ 2 . 

does not bring, in full accord with the scheme of Table 1, anything new, indeed. 
The general pattern indicated in Table 1 works and remains valid for all the matrix 
dimensions iV = AK, AK + 1, AK + 3 and AK + 3 with any auxiliary K = 1,2,.... 

On this basis, we are permitted to pay attention to any root Ct of our zero-order 
secular polynomial V (P) = (Cl — Q) <2(f2). As long as this root is an integer, it 
is extremely easy to study its small vicinity by the insertion of the ansatz Q = 
+ fi 2 u(fi 2 ) into the exact secular equation V(VL) = 0. The resulting new form of 
our secular equation reads 

p W) = ^ [^i(n) + . . . + f, 2K V K (Q)} . (29) 

On the basis of this formula we may consider a sequence ojj{ji 2 ) of approximate 
corrections initiated by u;_i(/z 2 ) = uo^ 2 {^ 2 ) — . . . — 0. The use of the abbreviations 
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Qj = Cl + fi 2 ujj(fi 2 ) re-interprets finally our secular equation (29) as an innovated and 
most important recurrent relation 

^ j " Qo(^-i) + - • ^ Qo(^-x + i) +/ ^ Qo(^-x) ' 3 ~ ' ' • • • • 

(30) 

At K = 1 and iV = 4 or iV = 5 this recipe returns us back to the current definition 
of the analytic continued fractions. Similarly, at N = 6 and iV = 7 we get their 
branched alterantive. At all the higher N > 8, our recurrences (30) offer just an 
immediate rational generalization of the latter two concepts. 



4 Interpretation and outlook 

4.1 Zero-order Schrodinger equation 

We have seen that once we switch to the rescaled energy variable E in eq. (11) and to 
the Hamiltonian (12), the zero-order diagonalization of H(0) becomes unexpectedly 
easy at any N. In spite of the manifest asymmetry of the Schrodinger QES equation 
in the limit £ — > oo, 





(N-l) 



1 





\ 



V 



(N — 1) 

1 



Po 
Pi 



\ 



PN-2 

\ Pn-i 



( 



E ■ 



Po 
Pi 



\ 



PN-2 

\ Pn-i ) 



all its eigenvalues remain strictly real, equal to integers and equidistant, 



(31) 



(Ei, E 2 , E 3 , . . . , En-i, E 



n 



'-N + 1, -N + 3, -N + 5, . . . , N - 3, N - 1) . (32) 



It is quite elementary to verify that also the the respective left and right eigenvectors 
of H(0) remain real. Up to their norm, all of them can be represented in terms of 
integers as well. Their components may be arranged in the rows and columns of the 
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following N by N square matrices P = P(N), 



P(l) 



P(2) 



V2 



p(3) = 7t 



i 



-1 



i 



1111 

2 0-2-4 

P(5) = ^6 6 "2 6 
' -2 2 -4 

-1 1-1 1 

All these matrices are asymmetric but idempotent, P 2 



1 

-1 

( \ 

3 
3 
1 

\ 



-1 -1 
-1 1 



1 1 N 
-1 -3 
3 



-1 



I. This implies that the 
Hamiltonian in our zero-order QES sextic Schrodinger equation H^p^ = fft^E^ 
need not be diagonalized at all. Indeed, as long as the iV— plets of the zero-order 
(lower-case) vectors ffi ^ are concatenated into the above-mentioned (upper-case) N 
by iV matrices P = P(°\ we may also collect all the pertaining eigenvalues in a 
diagonal matrix e^°\ In this arrangement, the unperturbed Hamiltonian is factorized, 
= Pe^P. As a consequence, the zero-order equation is an identity since, in 
our compactified notation, it reads Pe^PP = Pe^ and we know that P 2 = I. In 
the next section we show that and how similar notation may be used in all orders. 



4.2 Rayleigh-Schrodinger perturbation recipe revisited 

We have seen that 

At any finite value of the spatial dimension D we have seen in section 3 that the 
routine power-series ansatz of perturbation theory becomes applicable even though 
the unperturbed Hamiltonian H(fi) itself is non-diagonal. We may write 



H(\)=H<® + \HW, a 
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where the perturbation is an asymmetric one-diagonal matrix in our particular il- 
lustrative example. Even without the latter constraint we arrive at the textbook 
perturbative representation of our matrix Schrodinger equation, 

(H (°) + A ff «) • + A p« + . . . + A V X) + 0(X K+1 )) 

= (p(°) + . . . + AV^ + 0(\ K+1 )) ■ (EM + ... + \K E (K) + ( A AT+1)) . (33) 

The sets of the vectors for corrections pf \ j = 1, 2, . . . , N may be concatenated in 
the square matrices tyW. This enables us to re-write the first-order 0(A) part of eq. 
(33) in the particularly compact matrix form 

e (i) + P ^r(i) e (o) _ e (o) P $(i) = P # (i) p. (34) 
In the second order we get 

£ (2) + p ^(2) £ (0) _ £ (0)p ^(2) = p p(2) p + p #(1) ^(1) _ p ^(1) £ (1) (35) 

etc. This is a hierarchy of equations representing their source (33) order-by-order 
in A. Their new merit lies in their recurrent character. Their "old" or "input" 
data occur on the right-hand side of these equations, while the "new" or "unknown" 
quantities stand to the left. All the higher-order prescriptions have the same struc- 
ture. In all of them, the diagonal part of each equation (i.e., of (34) or (35) etc) 
determines the diagonal matrices containing energy corrections (i.e., or etc, 
respectively). Non-diagonal components of these matrix relations are to be under- 
stood as definitions of the eigenvectors, with an appropriate account of the well 
known normalization freedom which has been thoroughly discussed elsewhere [21]. 
All these relations just re-tell the story of our preceding section but after an appro- 
priate modification they may also be used for the evaluation of separate corrections 
in some more complicated QES systems [10]. 

In the conclusion, we may emphasize that the practical reliability of any pertur- 
bation prescription is, mostly, determined by the quality of the zero-order approxi- 
mation. In this sense, our present study offers also a broadening of their menu. In the 
light of our results one will be forced to make a more careful selection between the ES 
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and QES V^^d^l)- Indeed, in each of these respective extremes one encounters 
different difficulties. Still, the simpler ES-based choice of Vo(|f|) is predominantly 
preferred in practice. Almost without exceptions, such a decision is being made for 
one of the following three apparently good reasons. 

• Up to now, many QES-based constructions stayed within the mere lowest- 
order perturbation regime since the majority of the naive implementations of 
the Rayleigh-Schrodinger perturbation scheme becomes complicated in the high 
orders [8]. 

• A priori, any available set of the QES bound states ipo is, by definition, incom- 
plete. This makes the matrix form of the pertaining propagators non-diagonal 
and, hence, not too easy to use even in the lowest orders. 

• Last but not least, many ambitious QES models cannot be used as eligible for 
perturbations because they are almost frighteningly complicated even in zero 
order [6, 9]. 

Within our new perturbation prescription, all these shortcomings were at least par- 
tially weakened. At the same time, the appeal of the ES-based perturbation con- 
structions should at least slightly be re-evaluated as well: (i) the class of the available 
V"l- ES l(|f|) is really extremely narrow for many practical purposes; (ii) several impor- 
tant (e.g., double- well) models lead to a perturbation which is not "sufficiently small" 
in the sense of Kato [22]; (iii) one needs to go to the really very high perturbations 
in many cases of practical interest [23]. 

In such a setting we showed that the technical treatment of the QES-type zero- 
order approximations can be significantly simplified. This could help, say, in situa- 
tions where we have to deal with a phenomeno logical potential V(r) which cannot 
be well approximated by any available exactly solvable V^ £5 ](|f|) while there exists 
an exceedingly good approximation of V(r) by a partially solvable V^^d^l)- Un- 
der this assumption we may recommend the present construction as a guide to new 
applications of the old large— i expansion idea [24] , especially when a very high order 
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N characterizes the polynomial part of the QES wave function tp^{r), which would 
makes the simple-minded zero-order construction prohibitively complicated by itself. 

5 Summary 

The problem of construction of large QES multiplets has been addressed here within 
the framework of the so called large— £ expansion method. Its key idea is very 
popular and re-emerges whenever radial Schrodinger equation is considered at a 
large angular momentum I. Originally, this attracted our attention because in all 
these recipes there exists an obvious ambiguity in a free choice between alternative 
small parameters A = l/£ and \ s hifted = 1/(2 + P)- We should remind the reader 
that there exist in fact many alternative shifted-large— £ versions of the expansions 
which use this freedom in different ways [25] but, mostly, this parameter plays a 
certain not too essential variational role. In this context, the position of our present 
approach may be exceptional; in QES context the role of our "optimal" and unique 
P happened to be much more essential, influencing not only the rate of convergence 
but also, directly, the form of the perturbation series itself. 

Thus, although in most cases the transition to a shifted large— I expansion is being 
designed to extend its practical applicability, we imagined that the construction of 
the N— plets of quasi-exact sextic-oscillator bound states (r) does not admit any 
free choice of the shift P at all. On the contrary, our approach has been shown to 
prefer the unique, "optimal" definition of the shift P = P(N) which leads to one of 
the best available constructions of the QES spectrum. This is our main result. 

Our secondary motivation stemmed from the fact that whenever one moves be- 
yond the first few trivial multiplet dimensions N, the practical appeal of QES so- 
lutions fades away. In this sense we have verified once more that one of the most 
efficient strategies of solving "difficult" Schrodinger equation is provided by the stan- 
dard Rayleigh-Schrodinger perturbation series [22] and by many of its various prac- 
tical modifications and alternative implementations [4] . We succeeded in developing 
a prescription with several merits described throughout the text. 
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In a broader methodical setting we were guided by the idea that the quality 
of perturbative results depends on the two decisive factors, viz., on the number 
of the available perturbation corrections and on the rate of convergence of their 
series. In both these aspects our new approach passes the test very well. We have 
seen that in comparison with the standard approach summarized for our present 
purposes in section 2.2, we arrived at the more compact and, presumably, most 
easily generated representation of our QES energies. Moreover, in a sharp contrast 
to the manifest divergence of the former (though, admittedly, more universal) recipe 
[15], our algebraic approach represents the QES spectra by convergent infinite series. 

Our new implementation of some old ideas of perturbation theory looks promising 
and will definitely serve as a guide to the more complicated QES constructions in 
the future. For the sake of simplicity we have chosen just one of the simplest models 
where the sextic oscillator is not even complemented by a quartic force and admits 
merely a quadratic repulsion near the origin. This simplified many technicalities and 
strengthened our belief that any future transition to models with more parameters 
may and will proceed along the same lines. 
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Table 1. QES secular polynomials V{9) = V (n) + /i 2 Pi(fi) + . . . + fi 2K V K (n) 
(the trivial factor E ignored at odd N). 



dimension 


perturbation 


degree of Vj(VL) 


N 


in V(tt) 


7=0 


7 = 1 


j = 2 


7 = 3 


and 1 


absent (K — 0) 











2 and 3 


absent (if = 0) 


1 








4 and 5 


linear (if = 1) 


2 









6 and 7 


linear (if = 1) 


3 


1 






8 and 9 


quadratic (if = 2) 


4 


2 







10 and 11 


quadratic (if = 2) 


5 


3 


1 




12 and 13 


cubic (if = 3) 


6 


4 


2 
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